Large scale effective theory for cosmological bounces 
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An exactly solvable bounce model in loop quantum cosmology is identified which serves as a 
perturbative basis for realistic bounce scenarios. Its bouncing solutions are derived analytically, 
demonstrating why recent numerical simulations robustly led to smooth bounces under the as- 
sumption of semiclassicality. Several effects, easily included in a perturbative analysis, can however 
change this smoothness. The effective theory is not only applicable to such situations where numeri- 
cal techniques become highly involved but also allows one to discuss conceptual issues. For instance, 
consequences of the notoriously difficult physical inner product can be implemented at the effective 
level. This indicates that even physical predictions from full quantum gravity can be obtained from 
perturbative effective equations. 

PACS numbers: 98.80.Qc,04.60.Pp,98.80.Bp 



Our universe, extrapolated backwards in time, gets 
denser and denser, and eventually reaches such extremes 
that classical general relativity breaks down. This theo- 
retical limitation is to be solved by a quantum theory of 
gravity. Since fundamental quantum gravity implies dra- 
matic changes to our understanding of space and time, it 
can at best provide a non-intuitive, deeply quantum de- 
scription of the situation where even concepts of time 
break down. More intuitive pictures can be obtained 
semiclassically, where the only non-singular way a uni- 
verse can behave is as a bounce at small volume back 
to larger scales. Often, models with special, mostly ho- 
mogeneous, bouncing solutions in detailed descriptions 
are available, but their genericity or stability under per- 
turbations remains uncertain. Special properties are re- 
quired for such explicit descriptions which make bounces 
difficult to generalize in particular to inhomogeneities. 
In general, one can then only fall back to a fully quan- 
tum formulation. Nevertheless, under controlled pertur- 
bations bounce scenarios could survive and present cru- 
cial insights for cosmological scenarios. A self-consistent 
analysis of solutions and potential effects of perturbations 
is made possible by the effective description provided be- 
low in the framework of loop quantum gravity. 

We first consider the model of a free scalar field (f> in an 
isotropic cosmology, which will later be seen to play an 
important role for effective theory. The scale factor a of 
the metric then changes in proper time r according to the 
Fricdmann equation (a _1 da/dr) 2 = p (ignoring numeri- 
cal factors) where p = i(d0/dr) 2 is the energy density of 
4>. To quantize canonically we need to introduce canoni- 
cal variables given by the two pairs (c,p) with c = da/dr, 
\p\ = a 2 and (<fi,p<j>) with = a 3 d0/dr. The variable p 
can in general take both signs unlike the scale factor a 
because it also contains information on the orientation of 
space. Here, however, it is sufficient to assume positive 
p. The Friedmann equation then amounts to a condition 
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between the phase space variables. For a free scalar, 



P0 



/V2 = ±|cp| =: ±H 



(1) 



is a constant of motion, and we can view as 
internal time variable in which the dynamics of the sys- 
tem unfolds given by the Hamiltonian H. (The value 
p^ then plays the same role as constant energy values in 
classical mechanics and the sign in ([TJ selects the direc- 
tion of time.) A direct canonical quantization is not fully 
straightforward due to the absolute value in H — \cp\. 
Rather than using a quantization of this operator, one 
can quantize just cp and restrict states to the positive 
part of its spectrum. Since H is conserved, it suffices 
to pick an initial state which is a superposition of pos- 
itive eigenstates only. Moreover, we will be interested 
here in states which satisfy the semiclassicality condi- 
tion AH -c (H) for an initial state at large volume. 
Such states are sharply peaked around (H), which can 
easily be chosen positive. Thus, for such states we can 
work with a quantization of H = cp, dropping the abso- 
lute value. Then, using the symmetric operator ordering 
H = ^ (cp + pc) , we are lead to the Schrodinger equation 
dijj/dt = pdip/dp + Tpj) for a wave function ip(p, t). 

In quantum mechanics, one usually solves first for the 
wave function tp and then computes suitable expecta- 
tion values relevant for measurements. Here, it turns out 
to be faster and much more powerful to solve for such 
expectation values directly, using a form of generalized 
Heisenberg equations of motion. All information in the 
system is contained in expectation values (c), (p) and 
fluctuations or "rt-point functions," using the language 
of quantum field theory, 



G" 



((c-(c)r- a (p-m a ) s 



(2) 



\p\ = \\p\ 



for integer n > 2 and a = 0, . . . , n (using again symmetric 
ordering) [l[ . Their Heisenberg- type equations of motion 
are (£)' = ([c,H])/ih = (£), (pY - (%H])/iK = -(p) 
which is easily solved by (c)(t) = cqe*, (p)(t) — c^e -4 . 
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To proceed similarly for the G a ' n , we need one addi- 
tional ingredient since, e.g., G ' 2 = (c 2 ) — (c) 2 is not the 
expectation value of an operator and the commutator 
[(c) 2 , H] is not yet defined. We can simply do so refer- 
ing to the Leibniz rule to be satisfied for commutators, 
which gives [(c) 2 ,H}/ih = 2(c)([c,H])/ih = 2(c) 2 . By 
this procedure we have G°< 2 = ([c 2 - (c) 2 , H])/ih = 2G°< 2 , 
G 1 ' 2 = and G 22 = — 2G 2 ' 2 as well as similar equa- 
tions for G a ' n for n > 3. Also this is easily solved by 
G°< 2 (i) = c 3 e- 2 *, G 1 ' 2 ^) = c 4 and G 2 ' 2 (t) = c 5 e 2t . Fi- 
nally, fluctuations are subject to the uncertainty relation 

G°' 2 G 2 ^ 2 - \G^ 2 \ 2 > \h 2 \{[c,p})\ 2 (3) 

which implies C3C5 > S 2 /4 + c 2 . Our solutions repre- 
sent wave packets following exactly the classical trajec- 
tory with constant relative spread Ap/ (p) = V G 2 - 2 / (p). 

This model is thus exactly solvable in a very strong 
sense: we are not just able to obtain wave functions in 
closed form but the dynamics of all n-point functions de- 
couples and allows exact solutions. Such a behavior is 
realized in usual quantum mechanical systems only for 
very special solvable systems such as the harmonic os- 
cillator. Its similarity to our system is not accidental 
since both have quadratic, albeit different, Hamiltonians 
in canonical variables. Otherwise, non-linear coupling 
terms between n-point functions would result, not allow- 
ing analytical solutions. Such coupling terms can eas- 
ily be seen when, e.g., adding an anharmonic potential 
q 3 to the harmonic oscillator. Its expectation value can 
be written as (p 3 ) — (p) 3 + 6(p}G 2 ' 2 + 6G 3 ' 3 , demon- 
strating the coupling between (p) and G 2,2 . Calculating 
equations of motion as above shows that all of them are 
now coupled non-linearly. In intuitive terms, the cou- 
pling describes how a semiclassical wave packet spreads 
and deforms, back-reacting on the peak motion given by 
expectation values. 

A free scalar isotropic model thus serves as the analog 
of the harmonic oscillator (or of a free quantum field the- 
ory) in quantum cosmology. This observation is of cru- 
cial importance, since such free theories always provide 
the basis for perturbation theory in the form of effective 
equations [2,12]. Before discussing this we will first in- 
troduce changes implied by a loop quantization which we 
have not used so far. 

The classical system and the above quantization are 
singular because vanishing volume p = is approached 
arbitrarily closely by the solutions. Loop quantum grav- 
ity implies modifications to the dynamics [Jf, H| which lead 
to bounces provided that one starts (and, in this model, 
then remains) in the semiclassical regime 0. This has 
been studied in detail recently @, where the semiclassi- 
cality condition implied large values of the Hamiltonian 
H. A loop quantization leads to a modified Hamiltonian 



of the form /(p)sinc [3j, with a function 

f \p\(l + 0(£ 2 /p)) for|p|»£ 2 
IKP) \ \p\- n for|p|«£ 2 

depending on the Planck length £p and a positive n. In 
such an isotropic context, / is bounded from below by 
a positive number [3] in contrast to the classical func- 
tion p. If this is used in an "effective" Hamiltonian 
H = /(p)sinc, one can easily see how bounces occur 
because sine is bounded from above such that p can- 
not become arbitrarily small when it solves the equation 
H = f(p) sine for large constant H . Numerical solutions 
[6| to the resulting equation for wave functions are in fact 
surprisingly well described by the effective Hamiltonian 
psinc: a semiclassical peak follows the effective trajec- 
tory very precisely with negligible spread and deforma- 
tions of the wave packet in "time" <fi. This is certainly 
unexpected for a quantum system where in general wave 
packets spread quickly and move away from the classical 
trajectory. Our analytical solutions below will explain 
this behavior and prove the genericity of semiclassical 
bounces beyond what is possible in a numerical analy- 
sis. Unlike previous investigations, the systematic the- 
ory is well suited to a perturbation analysis which can be 
performed even if strong model-specific assumptions are 
relaxed. It thus provides the technology for developing 
realistic bounce scenarios. 

The expression p sin c, using f(p) ~ p for p»lp which 
as we will see is always the case for our solutions, is not 
quadratic in the variables (c,p) in contrast to the clas- 
sical Hamiltonian. This seems to remove the solvability 
observed before, as it generally happens when solvable 
Hamiltonians are modified. We are forced to use a dif- 
ferent Hamiltonian since the loop quantization does not 
allow one to act with c itself but only with the exponen- 
tials e ±IC appearing in the sine. To deal with the new 
Hamiltonian, we notice that H can be reformulated in 
new variables p and J := pe lc in such a way that it be- 
comes even linear. We then have H = — ^i(J — J) for a 
Hamiltonian in non-canonical variables (J,p) with Pois- 
son relations {p, J} — —iJ, {p, J} = iJ, {J, J} = 2ip. 
Thus, the Hamiltonian with these new variables forms a 
linear algebra isomorphic to sl(2,R). Classical equations 
of motion can easily be solved, also taking into account 
the reality condition J J = p 2 for our complex variables. 

If such a system remains linear after quantization, it 
allows a decoupling of the infinitely many quantum vari- 
ables to finitely many coupled linear equations just as 
Hamiltonians quadratic in canonical variables do. This 
would again make exact solutions possible. Quite sur- 
prisingly, the loop quantization H = —^i(J — J^) with 
basic commutation relations [p, J] — hJ, [p, ,V] = —KJ' 
and [J, jt] = —2hp — h 2 does remain a linear system. 

(These relations follow if the ordering J — pexp(ic) is un- 
derstood.) The classical algebra is simply modified by a 
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(trivial) central extension of charge h. In a quantum the- 
ory, there is an additional condition requiring that real 
variables are promoted to self-adjoint operators. Such 
conditions determine the physical inner product used to 
normalize wave functions. We do not use real variables, 
but the reality condition for J = pe lc implies JJ^ = p 2 
to be imposed (but note ,V J ^ p 2 in our ordering) . Im- 
posing the reality condition at the quantum level, which 
can easily be done for quantum variables as we will see 
below, corresponds to requiring that the physical inner 
product is used to normalize wave functions. 

Proceding as before, we derive generalized Heisenberg 
equations from commutators. This results in equations of 

motion (p)- = - 5«J) + (J)) and (J)' =-{p)-\h= (J) 

while (J) - (J) = 2iH is constant (where H = (H)). The 
general solution 

(p)(t) = !( Cie -* + C2 e*)-!?i (4) 
(J)(t) = i(cie"* - c 2 e*) + iH (5) 

differs from that obtained without a loop quantization 
and exhibits bouncing solutions (for c\C2 > 0) and sin- 
gular ones (ciC2 < 0). These solutions still need to be 
restricted by the reality condition J.P = p 2 . Taking ex- 
pectation values, this tells us 

\(j)\ 2 -{(p) + lK) 2 = \h 2 + G-G 2 > 2 (6) 

where we introduced the dispersion variables G 2 ' 2 — 
(p 2 )-(p) 2 &ndG=\{JJ^ + J^J)-\J\ 2 = (J.P)-\J\ 2 + 
hp+^h, using the commutation relations. As before, one 
can compute equations of motion for the dispersions and 
verify that the combination G — G 2 ' 2 =: e is constant in 
time. This implies a condition for the constants in our 
solutions p{t) and J(t): 

\(J)\ 2 - ((p) + ft/2) 2 = -4c lC2 + AH 2 = -e + h 2 /A . (7) 

Since e is a fluctuation term, it is small for semiclassical 
states, in which case cic 2 = H 2 + |e — j^h 2 is positive. 
In this case, ignoring small contributions from e and H, 
only bouncing solutions 

(p}(t) =Hcosh(t-8) , (J)(t) = -H(smh(t-S)+i) (8) 

with e s :— c\/H remain. Reality conditions of the effec- 
tive theory thus prove that all solutions which are semi- 
classical at at least one time are bounded away from the 
classical singularity at p = 0. (Note that this property 
of e <C H 2 is stronger than the condition AH <C H im- 
posed to take into account positivity of the Hamiltonian. 
There are thus physical solutions with large dispersions 
which reach the classical singularity, and the fact that 
states which are semiclassical at one time do not reach a 
singularity is non-trivial.) 

Similarly, we can explicitly determine the spread pa- 
rameters. We use the second order variables G 0,2 — 



(J 2 ) - (J) 2 , G 1 ' 2 = \{(pJ) + (Jp)) - (p)(J) and G 2,2 = 
(p 2 ) ~ (P) 2 - (The operator j\ although it appears 
in the Hamiltonian, is not included in quantum vari- 
ables here since we can use the commutation relations 
and the reality condition to express any expectation 
value containing Jt through quantum variables of the 
same or lower order not containing ,V .) By the gen- 
eral scheme we have equations of motion G ' 2 — —2G 12 , 
#,2 = _ 1go,2 - |G 2 ' 2 - \h 2 , G 2 - 2 - -2G 1 ' 2 . Thus, 
G ' 2 - G 2 ' 2 = 2c 3 is constant while 

G 2 - 2 = \{ Ci e- 2t + c 5 e 2t ) - ic 3 - ±h 2 
G — 2l c 4<5 — c 5 e ) . 

The constants c, are restricted again by the uncertainty 
relation which now, in non-canonical variables, reads 

G°< 2 G 2 < 2 - IG 1 ' 2 ! 2 > \h 2 \J\ 2 . (9) 

For t — > ±cx), this is easy to evaluate and satisfied only 
if C3, C4 and C5 are positive. Moreover, assuming H 3> h 
and near saturation before as well as after the bounce one 
has c A e~ 2S sa HH sa c 5 e 2S . With this, 

(Ap{t)) 2 = G 2 ' 2 (t) » HH cosh(2(t - 5)) (10) 

and not only the peak trajectory but also the spread of 
the wave packet is symmetric around t = 5. But only 
under the two stated conditions does the solution auto- 
matically become as coherent after the bounce as it was 
before (Fig. [J). 
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FIG. 1: Peak trajectory with surrounding area spread out by 
Ap around a bounce for H = 10ft. Dashed lines show the 
general behavior not assuming uniform spread. This figure is 
to be compared with the numerical solution in Fig. 2 of [5(- 

Our solutions in the given ordering are exact (20l | quan- 
tum solutions where (p)(t) is unaffected by spreading. 
This shows clearly why the naive effective formulation, 
obtained by simply replacing c 2 by sin 2 c in the classical 
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Hamiltonian, agrees so well with numerical simulations 
of a related free scalar model [f|. But the smoothness 
and symmetry of the bounce requires additional assump- 
tions which brings us to its robustness. As mentioned, 
the model is fully analogous to the harmonic oscillator 
as far as the decoupling of all its quantum variables is 
concerned. Using exclusively one solvable model to draw 
conclusions for the quantum nature of the big bang can 
then be as misleading as using the harmonic oscillator 
to determine the general behavior in quantum mechanics 
would be. Under any new ingredient to make the model 
more realistic coupling terms between all variables will 
arise. States then do spread and deform with compli- 
cated non-linear back-reaction on the expectation values. 
We thus need to discuss the expected more general be- 
havior, before coming to effective theory as a means to 
study stability and develop realistic scenarios. 

There are several examples for generalizations. First, 
while H 3> h is necessary for a semiclassical state at large 
volume of a homogeneous universe with large matter con- 
tent, several effects in inhomogeneous models remove this 
condition. Then, if is a spatial function and its local val- 
ues relevant for loop variables can be small even with a 
large matter content. (See for properties of inhomo- 
geneous states and operators defined in a regime around 
a homogeneous reference metric.) Moreover, with many 
more degrees of freedom, inhomogeneous situations allow 
the emergence of a semiclassical state at large volume, 
e.g. through decoherence, even though the state around 
the bounce can be highly non-semiclassical. This will 
complicate an analysis based only on the condition that 
the universe be semiclassical at one late time. 

But already in homogeneous models, several new ef- 
fects can arise if H is not large, as in those inhomoge- 
neous cases. The bounce scale will be smaller and quan- 
tum corrections to the p-dependence of the constraint 
through the function f(p) will appear. Since these cor- 
rections are higher powers in p, coupling terms between 
p and all quantum variables arise. At some point, due to 
the boundedness of f(p), no bounce occurs at all when 
H is small enough. (A similar removal of the bounce due 
to structural modifications to the equations, loosely re- 
lated to potential consequences of inhomogeneous states 
on small scales, has recently been explored in [13].) More- 
over, the behavior becomes more sensitive to the factor 
ordering of the constraint. Ordering all p to the right, 
for instance, leads to an additional contribution hp~ 1 J 
to the constraint. This is of the order h and would not 
change the above analysis much when H is large, but it 
does become relevant for smaller H. It implies coupling 
terms of the form —HG 1,2 /(p) 2 . In addition, there are 
corrections if model specifics are changed. Adding a cos- 
mological constant or positive curvature adds coupling 
terms to the Hamiltonian, but near the bounce they are 
not dominant. A potential for the scalar field, however, is 
more crucial because it leads to a time-dependent Hamil- 



tonian. 

In particular the latter change can have strong impli- 
cations depending on the form of the potential (see, e.g., 
the last example in [HI). The most important ingredi- 
ent to introduce is inhomogeneities and their evolution, 
at least as perturbations around a large isotropic back- 
ground, which is necessary both for a stability analysis 
and for cosmological predictions. Inhomogeneities add 
not only new degrees of freedom but also additional cou- 
pling terms between the variables. All this is difficult to 
analyze numerically, and generalized Hcisenberg equa- 
tions as used here can also become complicated to solve 
exactly. But all these changes are straightforwardly im- 
plemented in a general perturbation scheme around the 
solvable free scalar model identified here. The model be- 
havior is valid when only large semiclassical scales are in- 
volved and interactions are negligible, but it receives cor- 
rections when smaller scales are reached or interactions 
introduced. The resulting picture, accessible perturba- 
tively, can be quite different from the smooth bouncing 
solutions obtained for large H and isotropic states. 

This is not surprising because inhomogeneous situa- 
tions provide many degrees of freedom over which exci- 
tations can be distributed which happens in particular 
close to a classical singularity. The generic behavior is 
then different from an isotropic bounce, an observation 
which has also been made in the context of string the- 
121 ] . General singularity resolution can be achieved 



ory 

by mechanisms taking into account deep quantum be- 
havior which does not necessarily provide intuitive pic- 
tures IH, 14 , 15, lij] ■ While such a direct quantum anal- 
ysis is difficult in explicit terms, effective theory does 
make it feasible by perturbing around an exact model. 
Given a Hamiltonian H = Hq + Hi (p, J; t) with the solv- 
able one Hq and possibly time (i.e., scalar field) depen- 
dent coupling terms H\, we find perturbativc solutions 
(p) (t) = Po (t) + Pi (*) where po is a free solution as above 
and p\ is obtained by integrating 



Pi = -itr 1 {\p,H 1 ])\ ( ~ 



fill 



The right hand side is known in terms of zero order so- 
lutions. This procedure agrees with perturbation expan- 
sions in the interaction picture of quantum field theory. 
Indeed, effective equations derived by perturbation the- 
ory around our solvable model are in complete analogy 
with low energy effective actions in particle physics [J 0] • 
We thus provide the first systematic scheme for a well- 
defined evolution of perturbations through bounces. It 
allows us to approach the strong quantum regime and 
to determine if a semiclassical description would break 
down. One can, for instance, start with large H and 
then decrease its value to see corresponding changes in 
the wave function. (Although different in its nature, this 
parameter thus plays a role similar to others which have 
already been used in loop quantum cosmology [itI [l8| . 
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While the parameter should not be large for realistic in- 
homogeneous bounces, it provides a technical tool to sep- 
arate some quantum effects from others.) In perturbation 
theory one can then deal with quantum corrections aris- 
ing from several coupling terms which may be present. 

Effective theory is thus a powerful tool where direct 
numerical simulations become more involved. One can 
obtain directly quantities of interest such as expecta- 
tion values and fluctuations without taking the "detour" 
of wave functions. What is particularly interesting for 
full quantum gravity is that even conceptual issues can 
be addressed. The problem of time does not arise be- 
cause one is perturbing around a model which is explic- 
itly parameterized. Similarly, observables can be evalu- 
ated perturbatively. Even the notoriously difficult issue 
of the physical inner product, i.e. the issue of how to 
select normalizable wave functions, can be addressed at 
the effective level as already seen: we used reality con- 
ditions to select the physically viable effective solution 
which turned out to be bouncing. This corresponds to 
the correct normalizable wave function [f| . This solution 
is uniquely selected by referring only to quantum vari- 
ables side-stepping explicit states, an observation which 
is encouraging for effective developments in full quantum 
gravity. In effective equations, one can then implement 
the physical inner product order by order on effective so- 
lutions. In particular, this argument applied to general 
dynamical expressions provides the strongest indication 
so far that repulsive effects which have been observed 
throughout loop cosmological models are active even 
in a full setting. Although they are not strong enough 
to lead to general bounces at smaller scales of generic 
states, we have shown that perturbative regimes of inho- 
mogeneous bounces do exist. Whether or not this is suf- 
ficient for our universe can be tested in detailed analyses 
which are now being undertaken. Through effective the- 
ory, all issues relevant for physical predictions can thus 
be addressed. 
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